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Abstract 


This paper reports on the dynamic behavior of a 250 kW proton exchange membrane fuel cell power plant (PEM FCPP) and a 250 kW microturbine 
(MT) when operating in parallel. A load sharing control scheme is used to distribute the load equally between the PEM FCPP and the MT. For stand 
alone operation of a PEM FCPP, a set of batteries or ultracapacitors are needed in order to satisfy the power mismatch during transient periods. 
Using MT in parallel with the PEM FCPP helps in eliminating the need for storage devices. Models for the PEM FCPP and the MT with power, 
voltage and speed controls are used to determine the dynamic response of the system to a step change in the load. Simulation results indicate 
viability of parallel operation of the PEM FCPP and the MT. These results are obtained using MATLAB®, Simulink®, and SimPowerSystems®. 
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1. Introduction 


Due to strict power quality and system reliability require- 
ments, new approaches for power generation and transmission 
have evolved. One of these is the microgrid. Microgrids are low 
power generators clustered to supply residential and/or com- 
mercial loads. The primary advantage of the microgrid is the 
elimination of the cost and minimization of transmission sys- 
tem problems, since generators are built close to load centers. 
Also, the purchase and sale of energy is possible if the micro- 
grid system is grid connected. Fuel cells and microturbines are 
among the potential candidates that are suitable for supplying 
electric and thermal energy in microgrid system. 

Fuel cells are not only characterized by higher efficiency than 
conventional power plants, but they are also environmentally 
clean, have extremely low emission of oxides of nitrogen and 
sulfur and have very low noise. The components of the fuel cell 
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system include the fuel processing unit or the reformer, the fuel 
cell stack, and the power conditioning unit [1,2]. The reformer 
produces hydrogen to supply the fuel cell stack by processing 
any hydrocarbon fuel such as propane, methane, or methanol. 
Using the hydrogen supply and oxygen from the air, the fuel cell 
stack produces electricity and water through an electrochemical 
process. The output from the fuel cell is dc power. To provide 
power to a residential load, or to the electrical grid, a power con- 
ditioning unit is needed. Many models have been proposed to 
simulate the fuel cell in the literature [3—7]. Due to the low work- 
ing temperature (80-100 °C) and fast start up, proton exchange 
membrane (PEM) fuel cell power plants (FCPPs) are the best 
candidate for residential and commercial applications. Based on 
the model introduced in [3,4], a model for a 250 kW PEM FCPP 
is developed and used to study the dynamic behavior in response 
to a step change. 

Microturbines (MT) are small and simple gas turbines that 
have three main components: compressor, combustor, and tur- 
bine. The high pressure air from the compressor when mixed 
with the injected fuel forms a combustible mixture. The mix- 
ture is ignited in the combustor to produce hot gas flow, which 
is used to drive the turbine [8]. Two distinct types of MT are 
identified in the literature, the single- and split-shaft MT [9]. 
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Single-shaft MT is high speed in nature, where the electric gen- 
erator and turbine are mounted on the same shaft. Split-shaft 
MT uses a gearbox to connect the electric generator to the tur- 
bine shaft. Many models are used in the literature to describe 
the behavior of the MT [9-13]. In this paper, the model used 
in [10] has been modified by adding speed and voltage control 
loops. 

This study uses a 250 kW PEM FCPP connected in parallel 
with a 250kW MT to supply power to a common load. The 
paper focuses on the dynamic response of the system and the 
individual generators to a sharp increase and decrease in the 
load. 

The paper is organized as follows: Section 2 presents the 
configuration of the proposed system. In Section 3, a modified 
250 kW PEM FCPP model is introduced. Section 4 presents the 
MT model. Tests and results are discussed in Section 5 followed 
by the conclusion in Section 6. 


2. Parallel operation of PEM FCPP and MT 


Fig. 1 shows the connection of the MT and PEM FCPP in 
parallel. Many different configurations can be visualized based 
on the technology related to the MT and the PEM FCPP. In 
the proposed system, a synchronous generator (SG) is coupled 
to the MT shaft. Since the SG produces 660 V, a transformer 
is used to step down the voltage to 440 V to match the PEM 
FCPP voltage. The PEM FCPP produces a dc voltage of 440 V. 
A dc/dc converter with a gain of 1.5 is used to boost the dc 
voltage. A dc/ac inverter is used to convert the dc voltage to 
an ac voltage of 440 V. The SG and PEM FCPP ac voltages 
are controlled by controlling the field excitation and modulation 
index, respectively. The MT and the PEM FCPP are connected 
in parallel at the ac bus to supply a common load. Details of 
the MT and the PEM FCPP models and controls are given in 
subsequent sections. 

For parallel operation of the MT and PEM FCPP, a load shar- 
ing control is needed (Fig. 1). The load sharing control strategy 
as shown in Fig. 1 can be summarized as follows: 


e During steady-state conditions, the MT and the PEM FCPP 
share the load equally. 

e Due to the slower response of the PEM FCPP during transient 
periods, the MT compensates for the power mismatch until 
the PEM FCPP is able to supply its share of the load. 
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Fig. 1. Parallel operation of PEM FCPP and MT system. 
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Fig. 2. PEM FCPP system block diagram. 


3. PEM fuel cell model 
3.1. Model description 


In [3,4,14] a model of a PEM FCPP is introduced. In this 
paper the model has been modified to simulate a 250 kW PEM 
fuel cell. 

The model is based on simulating the relationship between 
the output voltage and partial pressures of hydrogen, oxygen, 
and water. The detailed model shown in Fig. 2 [3,4,14] includes 
reformer and power conditioning unit. 

The 250kW PEM FCPP model parameters are based on a 
440 V dc bus voltage with a stack current capacity of 94 A, 
and a cell voltage of 0.8 V. Based on the above figures, the 
PEM FCPP consists of six parallel stacks, each stack has 
550 cells in series. Using the indicated number of cells and 
stacks the 250kW PEM FCPP model parameters are given in 
Table 1. 


3.2. FCPP power control 


Power control scheme has been developed in [14] that can be 
summarized as follows: 


Pe = mVee Vs 


sin(ô) (1) 
where Pac is the ac power, m the modulation index, Veer the de 
voltage, V; the load voltage, X the external line reactance, and 6 
is the phase angle of the ac voltage mV. 


Table 1 

PEM FCPP model parameters 

Parameter Value 

Stack temperature 343K 

Faraday’s constant (F) 96,484,600 C kmol7! 
Universal gas constant (R) 8314.47 J (kmol K)~! 
No load cell voltage (Eo) 0.8 V 

Number of cells per stack (No) 550 

Number of stacks (Netack) 6 


K; constant = No/(4F) 
Utilization factor (U) 
Hydrogen valve constant (Ky, ) 
Water valve constant (Ky,0) 
Oxygen valve constant (Ko, ) 


1.4251 x 1076 kmol (s A)~! 


0.8 


4.22 x 1075 kmol (s atm)~! 
7.716 x 1076 kmol (s atm)” ! 
2.11 x 1075 kmol (s atm)~! 


Hydrogen time constant (ty, ) 3.37s 

Water time constant (tp 0) 18.418 s 
Oxygen time constant (To, ) 6.74s 
Reformer time constant (T1) 2s 
Reformer time constant (T2) 2s 
Reformer PI gain (C1) 0.25 
Conversion factor (CV) 2 

Activation voltage constant (B) 0.04777 A7! 
Activation voltage constant (C) 0.0136 V 
Internal resistance (Ri) 0.2778 Q 
External line reactance (X) 0.05 Q 

PI gain constants C2, C3 0.1, 10 
Voltage reference signal (V,) 1.0 p.u. 
Methane reference signal (Qmethref) 0.000015 kmol s~! 
Hydrogen-oxygen flow ratio (rH_0) 1.168 
Current delay time constant (Ty) 3s 


Assuming a lossless inverter: 


Pac = Pac = Veen! 


_ NstackNol 


qm 2FU 


(2) 


(3) 


where Pac is the dc power, / the stack current, gy, the input molar 
flow of hydrogen, Nstack the number of stacks, No the number of 
cells per stack, F the Faraday’s constant, and U is the utilization 
factor. 

From (1)-(3): 


Oe ees (4) 
sın SS R 7 
mV sNoNstack ae 


Assuming a small phase angle sin(ô) = ô, 


2FUX 


e (5) 
MV, NoNsua ? 


Eq. (5) describes the relationship between output voltage 
phase angle ô and hydrogen flow qy,. Eqs. (1) and (5) show 
that the active power as a function of the voltage phase angle ô 
can be controlled using the amount of hydrogen flow. 

Output voltage can be controlled by the modulation index m. 
The modulation index is controlled using a PI controller. The 
input to the PI controller is the error signal (difference between 
the ac terminal voltage Vac and reference voltage V,). 
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4. Microturbine model 
4.1. Model description 


As mentioned earlier, microturbines are classified in two 
categories, the single-shaft or high speed turbine and split- 
shaft or low speed turbine. In the single-shaft configuration, 
the compressor, turbine, and electric generator are mounted 
on the same shaft. The turbine speed is in the range of 
50,000-120,000 rpm. The frequency of the produced voltage 
is in the range of 1500-4000 Hz. To reduce the frequency to 
60 Hz, a cyclo-converter is used. In split-shaft microturbine, the 
electric generator is driven through a gearbox. The gearbox is 
used to reduce the speed to 3600rpm. Assuming a two-pole 
synchronous generator (SG), the frequency will be 60 Hz. In 
this case, no power electronic devices are needed for frequency 
conversion. 

In this paper, a split-shaft model is used to determine the 
dynamic behavior of a microturbine. In [10], the authors used 
the GAST model without speed control to simulate the split- 
shaft microturbine. The GAST model, which is developed by 
general electric (GE), is one of the most commonly used model to 
simulate gas turbines. In this paper, due to the use of aSG, a speed 
controller is developed and used with the GAST model. The 
model details and parameters are shown in Fig. 3 and Table 2, 
respectively. 


4.2. Simplified synchronous generator model 


A simplified model of a synchronous generator is given in 
[15]. A modified version of this model is presented in Fig. 4. 

The equations that are used to drive the synchronous machine 
model are explained in Appendix A. To simulate the syn- 
chronous machine, the authors used a predefined model existed 
in SimPowerSystems of the MATLAB software [16]. The model 
parameters are as shown in Table 3. 


4.3. MT speed control 


As shown in Fig. 3, speed control is achieved by comparing 
the rotor speed œw, with a reference speed wref. A PI controller 


O, 
> Duurbine 


AP 


Fig. 3. MT system block diagram with speed and power control. 
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Table 2 

MT model parameters 

Parameter Value 
Rated power (Prated-MT) 250 kW 
Real power reference (Pref) 1.0 p.u. 
Turbine damping (Dturbine) 0.03 
Fuel system lag time constant (T1) 10.0s 
Fuel system lag time constant (T2) O.1s 
Load limit time constant (T3) 3.0s 
Load limit (Lmax) 12 
Maximum value position (Vmax) 1.2 
Minimum value position (Vmin) -0.1 
Temperature control loop gain (Kr) 1.0 
Power control proportional gain (Kp) 0.1 
Power control integral gain (K;) 1.0 
Speed control proportional gain (Ks) 1000 
Speed control integral gain (Kx) 12.5 
Speed reference (ref) 1.0 p.u. 


Fig. 4. Synchronous generator block diagram with voltage control. 


is used to control the error signal (wref — wr). The PI controller 
output is then connected as a low value gate input. 


4.4. MT real power control 


Referring to Fig. 3, to control the mechanical output power 
from the microturbine, a PI controller is used. The input to the 
PI controller is the error signal (Prep — AP), where AP is the 
difference between the generated and the load power. 


Table 3 

Synchronous machine model parameters 

Parameter Value 
Rated power (Prated-SG) 250 kW 
Rated line to line voltage (Vratea) 660 V 
Frequency (F) 60 Hz 
Inertia constant (H) 0.822 s 
Damping factor (Kp) 33.63 p.u. 
Number of poles (P) 2 
Internal resistance (R) 0.02 p.u. 
Internal reactance (X) 0.3 p.u. 
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Fig. 5. Step change in load. 


4.5. MT voltage control 


As stated in Table 3, the SG output voltage is 660 V. To 
compare the SG voltage and the PEM FCPP voltage (440 V) 
a transformer is used to step down the SG voltage to 440 V 
at the load terminals. A PI controller is used to control the 
output voltage by controlling the excitation voltage of the syn- 
chronous generator. The input signal to the PI controller is the 
difference between the output voltage and reference voltage 
Vief. The parameters for the PI controller are Ky =0.005, and 
kK; =0.1. 


5. Test and results 


To evaluate the parallel operation of the PEM FCPP and 
the MT system, a step change in the load is used as illus- 
trated in Fig. 5. The simulation time is 100s, the initial load 
is 130kW, and increased to 260kW after 40s. After 70s, the 
load is decreased to 130 kW. The effect of changing the load on 
the output power of the PEM FCPP and the MT are shown in 
Figs. 6 and 7. 

As shown in Figs. 6 and 7, when the load increases the PEM 
FCPP is unable to supply its share instantaneously. In transient 
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Fig. 6. PEM FCPP output power. 
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Fig. 7. MT output power. 


condition such as a step-up in load demand, the microturbine 
is used to compensate the FCPP delay until the FCPP output 
power increases and matches half of the load. During the first 
40s, the system started with a load of 130 kW, initially, the MT 
provides a 130 kW while the FCPP output is zero. As the FCPP 
output power increases, the MT output power decreases. After, 
about 17s, the MT and the FCPP equally share the load. The 
delay in the FCPP output is mainly caused by the reformer, 
due to the slow gas processing response, as shown from the 
hydrogen flow and inverter phase angle (Figs. 8 and 9). The 
other reason of this delay is the thermal inertia, which plays an 
important role in the response time of the fuel cells. In a tran- 
sient condition, the temperatures of each stage of the reformer 
change slowly and continuously because of the high thermal 
inertia of the materials and the reactants. At t=40s, the load 
changes suddenly to 260 kW, the MT produces about 198 kW 
instantaneously after the load increases (Fig. 7). The MT output 
power deceases gradually as the PEM FCPP output increases. 
At time t=57s, the PEM FCPP and the MT share the load 
equally. 

During load reduction period, t>70s, the PEM FCPP is 
unable to decrease the output suddenly to 65 kW. In this case, 
the MT output at the load terminals decreases to —0.293 kW. 


Hydrogen flow (mol/s) 
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Fig. 8. Hydrogen flow in PEM FCPP. 
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Fig. 9. Inverter phase angle. 


This means that when ¢ is between 70 and 80s, the PEM 
FCPP is supplying the load as well as the losses in the trans- 
former and the transmission line between the SG and the 
load. 

Despite the fact that the MT response is much faster than that 
of the PEM FCPP, the mechanical power output from the turbine 
suffers a time delay as shown in Fig. 10. During such delay 
periods the power mismatch is compensated by the rotor inertia, 
which causes a momentary decrease/increase in the rotor speed 
as illustrated in Fig. 11. During the low/high rotor speed periods, 
the speed controller and the power controller increase/decrease 
the input power to the turbine. When the mechanical input power 
matches the electrical power output, the speed controller brings 
the rotor speed back to synchronous value (Fig. 11). As indicated 
in Fig. 12, the summation of the output power from the PEM 
FCPP and the MT matches the step change in load as given in 
Fig. 5. 

The change of the ac voltage of the PEM FCPP is insignif- 
icant (Fig. 13). This is mainly due to the effectiveness of the 
PEM FCPP voltage controller as shown in Fig. 2. Fig. 14 shows 
the change of fuel cell stack dc voltage due to load changes. 
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Fig. 10. Turbine mechanical power. 
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Fig. 11. Synchronous generator rotor speed. 
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Fig. 12. Total output power from the PEM FCPP and the MT. 


In order to produce a constant ac voltage, the modulation index 
changes with inverse proportionality with respect to the dc volt- 
age (Fig. 15). 

During load changes, the MT voltage is experiencing a short 
period of oscillation as illustrated in Fig. 16. As a result of the 
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Fig. 13. PEM FCPP ac voltage. 
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Fig. 15. Inverter modulation index. 


voltage oscillations and load changes, the voltage controller 
changes the excitation level to keep the output voltage con- 
stant at 660 V. Fig. 17 shows the changes in the excitation 
voltage. 
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Fig. 16. MT voltage. 
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Fig. 17. SG excitation voltage. 


6. Conclusion 


Dynamic plots reflecting the parallel operation of PEM FCPP 
and MT are evaluated. Dynamic models of the 250kW PEM 
FCPP and MT with speed and power control are developed. 
The step load change test results indicate that the proposed 
system configuration is able to accommodate drastic changes 
in the load. In stand alone operation of the PEM FCPP, a 
set of batteries or ultracapacitors are essential in order to sat- 
isfy the power mismatch during the transient period. Using the 
MT in parallel with the PEM FCPP can help in eliminating 
the need for storage devices. Although the MT response has 
a better load following capability, the mechanical power out- 
put experiences a momentary delay. During this delay period, 
rotor speed decreases/increases as part of the rotor kinetic 
energy is used to compensate the power mismatch. By the 
time the MT system reaches a new equilibrium point, the 
speed controller brings the rotor speed back to synchronous 
speed. To generalize the fact that the MT can be used to 
preclude the storage devices in FCPP, more simulation for 
different situations and different control scenarios may be 
used. 
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Appendix A 


Synchronous machine model as explained in [15] is con- 
structed based on synchronous machine swing equation. Swing 
equation for power system dynamics can be expressed in many 
forms. Given @ in electrical deg s71, swing equation is expressed 
as: 


(A.1) 


Linearizing the swing Eq. (A.1) results in an equation as 


shown in (A.2). 
tjAws = AT, — AT, (A.2) 


The basic equations for the simplified linear synchronous 
generator model consist of three equations, i.e., 


; K3 K3K4 

AE, = ——>—— AEpp - ———— (A.3) 
q 1+ K3t 98 1+ K3t 98 

AT, = Kj Aô + K2AE, (A.4) 

AV, = KsA6b + KoAE,. (A.5) 


The constants K1, K2, K3, Ka, K5, and Ke are dependent upon 
the network parameters, the quiescent operating conditions, and 
the infinite bus voltage. 

The simplified synchronous generation model as presented in 
Fig. 3 is obtained by combining the effect of the damping torque 
with voltage control, and Eqs. (A.2)-(A.5). The variables are 
defined as follows: 


Tm mechanical torque 

Ta accelerating torque 

Te electromagnetic torque 

Tj time constant (t; =2H/a,) 

Tho d-axis transient open circuit time constant 

ô rotor angle 

Vi terminal voltage 

Ey rms of the peak stator voltage 

Erp stator EMF 

Or rotor angular speed 

Kı change in the electrical torque for small change in rotor 
angle at constant d-axis flux linkage 

K2 change in the electrical torque for small change in the 
d-axis flux linkage at constant rotor angle 

K3 impedance factor 

K4 constant related to demagnetizing effect 

Ks change in the terminal voltage V; for a small change in 
rotor angle at constant d-axis flux linkage 

K6 change in the terminal voltage V; for a small change in 
the d-axis flux linkage at constant rotor angle 

Kp damping constant 

H inertia constant 
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